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Abstract 

Using data obtained through simultaneous measurements of {p,p') spin- 
transfer observables and {p,p'^) coincident spin observables, we have made 
a model-independent determination of the complete scattering amplitude for 
the 15.11 MeV, 1^, T=l state in ^^C at an incident proton energy of 200 MeV, 
for four proton scattering angles ranging from ^c.m. = 5.5° to 16.5°. At each 
angle, 16 different observables were determined, whereas only 11 independent 
quantities are required to specify the transition amplitude for this state. It 
had been shown previously that the set of observables measured span the 
allowed space; hence the system is over deter mined, which allowed us to ex- 
tract, in a model-independent fashion, each of the individual spin-operator 
amplitudes that characterize the reaction. Additional insight into the phys- 
ical mechanisms that drive this transition is obtained by mapping out the 
momentum-transfer dependence of these amplitudes. We also compare the 
magnitudes and phases determined for each of the spin-operator amplitudes 
to the predictions of calculations performed in both relativistic and nonrela- 
tivistic frameworks, and discuss the physics content of these comparisons. 

25.40.Ep, 24.70.-hs, 25.90.-hk 
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I. INTRODUCTION 



At intermediate energies of ~150-500 MeV, hadron-induced nuclear reactions have served 
as a rich source of information, both to further our understanding of nuclear structure, and 
to illustrate how the elementary nucleon-nucleon {NN) interaction may be modified in 
the nuclear medium This is particularly true of spin observables, which often contain 
contributions from interferences between different pieces of the full scattering amplitude. 
These interferences, inaccessible through measurements of differential cross sections alone, 
can thus provide information about not only the magnitudes, but also the relative phases 
of individual terms in the transition amplitude. Ideally, one would like to extract precise 
values for each of these (complex) terms in a manner which does not rely on specific model 
assumptions. This information, in turn, can then provide the most stringent tests of a given 
theoretical prediction for the hadronic process of interest. 

To obtain deeper insight into the nucleon-nucleus (NA) interaction, several inelastic 
transitions have been identified as being particularly amenable to experimental investiga- 
tion. The strong, isovector 1"^ state in ^^C at an excitation energy of 15.11 MeV has been 
extensively studied, using a variety of probes operating under a broad range of kinematic 
conditions. Due to its l"*" nature, hadronic excitation of this state represents an "unnatural 
parity" transition; consequently, this process should be particularly sensitive to the spin- 
dependent parts of the NA interaction, and in fact has long served as a critical test for our 
understanding of the AS = AT = 1 component of the NN effective interaction 0. As an 
example, the spin observable P — Ay, i.e., the difference between the polarization induced in 
the outgoing proton and the scattering yield asymmetry that results from use of a polarized 
incident beam, has been investigated both experimentally and theoretically 0. This 
quantity, which in a non-relativistic framework contains only interference terms between 
competing pieces of the transition amplitude, has been shown to be sensitive to the cou- 
pling of the nucleon spin to the bound nucleon current ^ . In later relativistic treatments of 
proton-nucleus scattering [^|^ , it became clear that these same nuclear current terms appear 
more naturally in a relativistic formalism, and are produced through the linear couplings 
between the upper and lower components of the bound nucleon wave function [^|^ . In either 
description, P — Ay is dependent on the momentum of the nucleons inside the nucleus, and is 
therefore sensitive to the off-shell behavior of these nucleons, and the non-local or exchange 
nature of the interaction. 

Thus, detailed comparisons of measured spin observable data to various predictions for 
specific nuclear transitions has traditionally been our best means of constraining theoretical 
models. In recent years, though, such comparisons have raised a number of concerns, at 
both the experimental and theoretical ends. For example, it has been shown |^,|^ that in a 
direct-only plane- wave impulse approximation (PWIA), certain combinations of {p,p') spin- 
transfer observables are directly related to individual terms in the effective A^A^ interaction, 
weighted by corresponding nuclear response functions. When physical processes such as 
distortion of the projectile waves and knock-on exchange are included, however, it becomes 
less clear what these combinations of {p, p') observables represent physically, thereby blurring 
any simple interpretation of these quantities. On the experimental side, counting arguments 
show that even "complete" sets of {p,p') spin-transfer measurements provide (at most) 
eight independent pieces of information, while a transition of the general form 0"*" J'^ 
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requires knowledge of (8J + 3) independent quantities in order to fully specify the scattering 
amplitude. 

To reduce this ambiguity, it is clearly desirable to establish a more direct connection 
between theory and experiment, such as would be obtained from a model- independent ex- 
traction of individual elements of the scattering amplitude. In so doing, one bypasses the 
'intermediate' role played by spin observables, which represent non-trivial combinations of 
many of these elements. Carrying out such a program, though, requires access to informa- 
tion beyond that provided by {p,p') spin-transfer measurements alone. In particular, one 
must probe the polarization state of the 'other particle' involved in the reaction: the recoil 
nucleus. This can be achieved through study of the angular correlation in the final state 
between the scattered (out-going) proton and the particle (s) emitted in the decay of the 
excited nucleus, as in reactions of the type (p,p'7). More specifically, for proton excitation 
of a O"*" — i> l"*" transition, followed by electromagnetic decay back to the ground state, it has 
recently been proven formally [|T3,0 that certain {p,p''j) measurements, when combined 
with complete sets of {p,p') observables, allow for a complete description of the transition 
amplitude. Note that in this case a total of 8J + 3 = 11 independent quantities must be 
determined. 

In this paper, we report the first model-independent determination of the complete tran- 
sition amplitude for any nuclear final state with J 7^ 0. This analysis is based on a detailed 
set of measurements ||12| carried out at the Indiana University Cyclotron Facility (lUCF), 
in which both {p,p') spin-transfer and {p,p''y) coincident observables were determined si- 
multaneously for excitation of the isovector 1+ state at 15.11 MeV in ^^C. Data were taken 
at four proton scattering angles (6'c.m. = 5.5°, 8.8°, 12.1°, and 16.5°) at an incident pro- 
ton beam energy of 200 MeV. As will be shown below, this data set allows for extraction 
of each individual (complex) spin-dependent term in the transition amplitude in a model- 
independent manner. Moreover, because this procedure has been applied at a number of 
scattering angles, the momentum-transfer (g) dependence of each spin-operator amplitude 
has been mapped out, allowing for a clearer interpretation of the physical mechanisms that 
drive this transition. 

This paper will be organized in the following way. In Sec. II we present the formalism 
adopted (from Ref. [T^) for analysis of the data. In Sec. Ill, we present the observables mea- 
sured in the lUCF experiment [|I2|, and express each of these in terms of real and imaginary 
parts of the individual spin-operator amplitudes. We describe in detail the procedure used 
to extract these amplitudes from the observables in Sec. IV, and present the results of this 
procedure in Sec. V. Section VI provides a detailed discussion of the momentum-transfer 
dependence determined for each amplitude (both magnitude and phase), and compares our 
results to predictions of calculations performed in both relativistic and nonrelativistic for- 
malisms. Where possible, we discuss the specific physics issues to which each amplitude is 
sensitive. Our most significant results and conclusions are summarized in Sec. VII. 

II. FORMALISM 

The formalism adopted here follows that of Piekarewicz et al. [|T^ , and has been presented 
in detail in Ref. [|12|; only the most important results are given here for completeness. We 
work in the orthogonal coordinate system defined by 
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n = pxp'; K = p + p'; q = nxK, 



(1) 



where p (p') is the incident (outgoing) proton momentum, n is directed normal to the 
scattering plane, K is along the direction of the average proton momentum, and (neglecting 
the reaction Q-value) q points in the direction of momentum transfer p'-p. In this frame, 
the most general form of the scattering amplitude allowed by angular momentum and parity 



conservation for a ^ 1 transition can be written in the form [10 



fP(p,p') = A„o(S-n) + A„„(S-n)(a-n) + A;,k(S ■ K)(a • K) 

+ A^,(S • K)((T ■ q) + A,^(S ■ q)(a • K) + A,,(S • q)(a ■ q) , (2) 



where 



|1+M)(0+| (3) 



is the polarization operator for promoting a O"*" state to a state with J'^=l+ and magnetic 
substate M, while a are the Pauli spin operators for the projectile. In Eq. 0, the Aj^'s 
are scalar functions of energy and momentum transfer, where the subscripts i{= n, K, q) 
and /i(= 0,n,K,q) indicate the polarization components of the recoil l"*" nucleus (S) and 
the scattered proton (a, with ao = 1), respectively. Because there are six allowed complex 
amplitudes, 11 pieces of information (after eliminating one overall phase) are required to 
specify the scattering amplitude for this transition. 

If the final polarization of the nuclear state is undetected, one can sum incoherently over 
the S index. In this case, the relevant spin observables can be expressed as 

_ ^A. A* X.. 



d(7 1 

Dap = AfiA*^6ij-Tr{(7a<7^(7i3au} , (4) 



dilr, 

where a,(3 = 0,n,K,q. Using the scattering amplitude of Eq. ^ and carrying through the 
Pauli algebra, it is easily shown that only eight of the 16 possible singles {p,p') spin ob- 
servables (the Da/3) are nonzero. Because 11 independent quantities appear in T^(p,p'), 
one sees that singles measurements alone can not uniquely define this amplitude; or, put 
another way, there is information contained in TP{p,p') which is not accessible via {p,p') 
observables. In particular, the presence in Eq. | of the Kronecker 5ij makes it impossi- 
ble to determine the relative phase between any two spin operator amplitudes (Aj^'s) that 
correspond to orthogonal orientations of the recoil nuclear polarization. 

If one makes the assumption that the {p,p''y) reaction is strictly a two-step process, 
then the total transition amplitude (excitation plus decay) can be written as the product 
of the strong and electromagnetic amplitudes. In this case, all of the coincident {p,p''y) 
observables may be written in terms of just the singles {p,p') spin operator amplitudes, ^j^, 
and the 7-ray branching ratio to the ground state R fl^. In complete analogy with the 
singles observables (Eq. the spin-dependent coincident observables for this transition can 
be written in the form 
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^'"^ (k) = |?E^^.^y.(k) 



dVtpdVt^ Svr f^^ 

wo Jo (k)-Pa/3(k) = ^ Ai^A*^tij(k) ^Tr{a„(T^a^(T^} , (5) 

where k is the momentum direction of the emitted photon, and the photon polarization 
tensor tij is given by 

= 6,j - (k ■ ei)(k ■ e^) , (6) 

with §4 a unit vector lying along one of the (n, K, q) coordinates defined in Eq. |l|. It is 
important to emphasize that the coincident observables defined in Eq. ^ are written in terms 
of the appear in the definitions of the singles observables, Eq. ^. The crucial 

difference between these two sets of equations is the presence of tij(k) in the definition of the 
coincident observables, replacing the Sij for the singles observables. Thus, certain [pjp'j) 
spin observables will be sensitive to the relative phases between amplitudes for different 
recoil nuclear polarizations, provided that the emitted 7-ray has momentum components 
along both ej and ej. Thus, the specific information about the transition amplitude which 
is lost for the singles observables can be accessed through coincident [p, p''~f) measurements. 



III. SUMMARY OF OBSERVABLES 



In this paper, we present detailed analysis of the data from a previous experiment |]12 
In that work, some spin observables were measured using an incident proton beam whose 
polarization vector was normal to the (horizontal) scattering plane, while others required 
use of a beam in which the polarization vector had been precessed to lie in the scattering 
plane. Due to technical problems encountered during data acquisition [T^, reliable values of 



the differential cross section, da/dilp, could not be extracted from the experimental yields. 
While this problem had no effect on the determination of spin-dependent observables, it 
was necessary to use previously measured values of da/dflp for this transition [|l^ to set the 



scale for the overall magnitudes of the y4j„ amplitudes. In this paper, these amplitudes will 



be presented in units of y/ib/sr. 

The singles observables obtained in Ref. [l^ include the three normal-component spin- 
transfer coefficients Ay, P, and Dn'n-, and two linear combinations of the in-plane spin 
transfer coefficients, 

Dx = Dl'l sin a + Ds'i^osa , 

where here a is the spin precession angle about n experienced by the scattered protons in 
passing through the magnetic spectrometer used in the experiment. Through use of Eq. ^, 
all of these measured (p, p') observables can be expressed in terms of the Ai^ amplitudes as 
follows: 
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= 2[^{A,^,A*^^)-'^{AkkA*k, + A,kA*^^)], 



(Kip 



and 



-^Dx = cos9pK[sm{a-9c.m.+dpK)-^DKK + cos{a-9c.m. + dpK)-^DgK] 
= sm9pK[sm{a-9c.m.+9pK)-j7z-DKK + cos{a-9c,^, + 9pK)-jT-DgK\ 



dQp " ' "'"dQp '''' ^ ^-"^ "''''dQp 

+ cos 9pK[ sin(Q; — ^c.m. ~l" ^pi<r) ^^ D]^q -\- COs(tt — ^c.m. ~\~ 9pK) -^gg] ' l"^) 

where 6*0.111. is the proton center-of-mass scattering angle, 9pK is the angle between the incident 
beam momentum p and the average momentum K (see Eq. ||), and a is the spin precession 
angle defined above. The Dap spin transfer coefficients that appear in the previous equation 
can also be expressed in terms of the spin operator amplitudes as: 



dVtp 



—Dk, = -2[%A^oAl^)-^{AKKA*Kg + A,KAlg)], 



dflp 



'^'^ DgK = 2[S(A„oA:J + mKKA*j,g + AgKA*J] , 



dflp 



^^-Dgg — |yl„oP — l^nnl^ — l^ii'if 1^ + l^i^gl^ — lAgE'P + l^ggl^ • (10) 

For the coincident [pjp'j) measurements, the incident beam polarization was either 
pointed along the normal n to the scattering plane, or was rotated to lie in the scatter- 
ing plane. During vertical polarization running, three photon detectors were positioned in 
the scattering plane on beam right, and one was placed directly above the target [|l^]. If 
a sum over the two beam spin states was performed (effectively producing the results that 
would be obtained with an unpolarized beam), the coincident yields from the three in-plane 
detectors could be used to extract the in-plane double-differential cross section, which takes 



the form 10 



r2 

^ ^ = A{9p) + B{9p)cos29^ + C{9p)sm29^ , (11) 



3R dVt^dVtp 
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where R is again the 7-ray branching ratio to the ground state, 6^ is the photon angle in the 
scattering plane, measured with respect to the q direction, and A, B, and C are unknown 
functions of the proton scattering angle. These can be written in terms of amplitudes as: 



A = \Ano\^ + \A^n\^ + ^[\Akk\^ + \AKg\^ + \AgK\^ + \Agg\ 
1 



B = -[\Akk\'^ + lAxql'^ — \AgK\'^ — lAqql"^] , 

C = -^AkkAIj, + AKgA*J . (12) 

We point out that in this work our definitions of A, B, and C do not include the branching 
ratio normalization factor 87r/(3i?), and are therefore different from the definitions given in 
Ref. |T^. A more important comment is that, in order to eliminate any sensitivity to detec- 
tion efficiency in the photon detectors, only the ratios of the symmetric and antisymmetric 
pieces of the in-plane coincident cross section, B/A and C/A, were used in our determi- 
nation of the transition amphtude. These quantities have the obvious advantage of being 
independent of any overall normalization error. 

We now turn to the spin-dependent coincident observables which, at a given proton 
scattering angle 6p, will be functions of the angle 6^ at which the photon is emitted. If the 
photon is emitted in the scattering plane, then the normal-component coincident analyzing 



power, scaled by the coincident cross section, can be cast in a form similar to that of Eq. |TT 
i.e., 

"^"^ ie^)Ay{e^) = eA{0p)+eB{9p)cos29^ + ec{9p)sm29^ , (13) 



3R dVt^dVtp 



where 



= 2g?(A„oA:j + ^{AKKA*Kq + AgKA*qq) , 

= ^{AKKA*^q - AqKA^q) , 
— _(^( A A* _ „ A* \ 



(14) 



tB — ^\n.KK^Kq~ ^qK^qq) ■ 
eC = -^{AKxAlq - AxqA^j^ 



and d'^a/dQjdQp is given by Eqs. |TT] and |T2|. Note that with three values measured for 
the normal-component spin asymmetry (corresponding to the three angles of the photon 
detectors), we have three expressions from which the coefficients e^, e^, and ec can be 
determined. Although this inversion is possible in principle, it is more efficient to use the 
measured asymmetries directly in our fitting procedure to determine the independent Aj^'s, 
which we describe in the next section. 

With a normally-polarized beam, the photon detector directly above the target could 
measure another piece of the normal- component coincident analyzing power, given by 

n)A,(n) = 2^{AKKA*j,q + AqKAlg) , (15) 



dQ dQ ~ ^^y^KK^Kq^ ^qK^qq 



where 

(n) = \Akk? + \AKq? + \AqK\' + \Agg\' . (16) 



dVt^dVtp 
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It can be shown algebraically that the ratio of Eqs. ^ and can be expressed in terms 
of the normal-component singles observables in the form 

^.(fi) = - }f~rf'\ ■ (17) 

— l>n'n) 

In Ref. this equation was used to demonstrate that the independently measured {p, p'7) 
coincident and {p,p') singles observables were internally consistent. 

Finally, the four remaining (p, p'7) observables were extracted from the spin asymme- 
tries measured when the incident beam polarization had been oriented to lie in the scattering 
plane. These are components of the longitudinal and sideways analyzing powers, DoL(k) 
and Dos'(k), respectively. These quantities, which vanish identically for singles measure- 
ments, are related to the corresponding center-of-mass asymmetries Doir(k) and Dog(k.) via 
a rotation in the reaction plane through the angle Opx'- 



fDoL{k)\ ^ /cosV -sinVU^ox(k)^ 
\^Dos (k) / \sm 9pK cos 9pK / V ^og (k) y 



(18) 



with Dqk{}^) and -Dog(k) defined by 



a 



3R dVtpdVt^ 



-(k)Z}ox(k) = 2[g?(A„oA^^)-53(A„„Ay]W(k) 

+ 2[3f?(A„o^:i,) - %A^nAl^)]tr,,{y) , 

^^^^(k)/^o.(k) = 2[5J(A„oA^J + 53(A„„A;,^)]W(k) 

+ 2[3ft(A„oA;^) + 53(A„,„A;^)]t„,(k) . (19) 

The four coefficients of t„i^(k) and t„g(k) that appear in the above expressions (i.e., these 
four combinations of the A^^'s) can be viewed as independent observables, and were deter- 
mined in Ref. |12| by a fit to the dual sinusoidal dependence of the measured asymmetries on 



both the outgoing photon direction k and the orientation of the incident proton polarization 
in the scattering plane. Contained in these four observables are clear sensitivities to the 
relative phases between terms in the transition amplitude corresponding to proton and re- 
coil nuclear polarization projections that are normal to, and oriented in, the reaction plane; 
relative phases which, by virtue of Eq. |[ are inaccessible via singles (p, p') measurements. 

IV. MINIMIZATION PROCEDURES 

In the set of combined {p,p') and {p,p''y) observables discussed above, we have a total of 
16 measured quantities at each proton scattering angle, which can be used to determine the 
11 independent quantities required to specify the complete scattering amplitude. Although 



it has been shown formally |11] that certain combinations of these observables can provide an 
analytic solution of this problem (via matrix inversion), this method of analysis has several 
disadvantages in practice. Because only particular linear combinations of the measured 
quantities are used, some statistical information is invariably lost in this method. Of more 
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concern is the loss of any statistical gauge of the quality or internal consistency of the 
amplitude extraction procedure; because one will always obtain an 'answer,' the assignment 
of errors to the Ai^ amplitudes becomes somewhat ambiguous. 

In this work, all of the observables measured at each scattering angle are used as input 
to a single minimization procedure. This method not only makes use of the full statistical 
information contained in the data set, but also minimizes any bias that could be introduced 
through data manipulation (e.g., forming various linear combinations of the data) prior 
to the actual fitting. Thus, our approach can be viewed conceptually as the following: 
we seek values for the six complex amplitudes Ai^ contained in Eq. ^ that, when used to 
form the 16 combinations provided in Section II, minimize the differences between these 
combinations and the measured values of the corresponding observables. Details of the 
minimization procedure we followed will be provided in the next few paragraphs, but we 
first point out a few subtle features unique to this problem. Because the amplitudes are 
complex, one can seek values for either their magnitudes and phases, or for their real and 
imaginary components. In this work, we carried out independent fits to both parameter 
sets, and obtained (as one would hope) equivalent results. However, it became clear that 
when attempting to resolve various discrete ambiguities in the fitting results, or invoking 
arguments of 'smoothness' (in momentum transfer), use of magnitudes and phases as the 
fitting parameters was favored. It was also necessary to hold one phase fixed during the 
fitting, as the observables are insensitive to any uniform shift in all the phases. Because 
\Ann\ was consistently one of the largest magnitudes over the entire angular range studied, 
we chose to define Ann to be real and positive, and thus determined (in effect) the phase of 
every other Ai^ relative to that of Ann- 

A final choice which required careful thought was the selection of appropriate starting 
values for the parameters to be fit. To avoid any bias, we noted that the form of the singles 
cross section, da/dflp (Eq. H), reveals that the maximum allowed value for the magnitude 
of any of the Aj^'s is constrained by 



Thus, the allowed parameter space for the magnitudes must be restricted to this range. In 
our fitting, we therefore assigned starting values for each by using a linear random 
number generator to select a number between zero and this maximum. Similarly, starting 
values for each phase were chosen at random over the allowed range 



By holding 0„„ at zero, then searching the 11-dimensional parameter space for minima in 
X^, we were able to determine the values for the that were most consistent with our 
entire 16-observable data set. 

We now describe the actual fitting procedure in greater detail. We begin by defining 



Ai,,\ ^ 



da 



(20) 




(21) 



16 



(22) 



j,k=i 
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where Fj are the measured values for the 16 observables, and fj are the expressions given in 
the previous section for these same observables in terms of the Aj^'s. In this equation, Wjk 
is the weight matrix, given by 



jk 



(23) 



Here e is the full error matrix associated with the set of observables. In the absence of 
any correlations among the observables, i.e., if each of the observables had been determined 
independently, Wjk would be diagonal, with each element equal to the inverse of the square 
of the error assigned to each observable. 



W. 



n 



5Fj 



(24) 



In this work, however, Wjk has been generalized to include known correlations (off-diagonal 
elements) between specific observables, given their method of determination 0]. For exam- 
ple, the values obtained for the in-plane coincident cross-section coefficients, B/A and C /A, 
were deduced from the same data set by fitting sinusoidal functions to the measured photon 
yields. Thus, the resulting coefficients of the fit are highly correlated. 



The minimization procedure we employed uses a combination of algorithms |TJ] to locate 
the minima of an arbitrary, nonlinear function in an arbitrarily large parameter space. The 
function (Eq. ^) can be linearized around a minimum value Xo 
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X 



Xo + E 



,=1 dA, 



(25) 



In this and following equations, we have simplified our notation by denoting all Ai^ with a 
single index, A^, representing either a magnitude or phase. In the above equation, Aq is the 

Xo- Thus, our minimization 



set of A^ values which minimize 
condition can be defined by 



i.e., their values at x^ 



dx^ 



dx^ 



dA,^ 
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d'x' 



+ E 



SA„ 



. 



(26) 



Once initial values for the A^'s are chosen, the parameter search begins by solving for optimal 
changes in the parameters, SA^, then generating new values for these parameters via 



(27) 



from which point the parameter search can begin again. This procedure is repeated until a 
minimum in x^ is found. 

The minimization condition can be written as a matrix equation 



a 5 J 



/5 



(28) 



or 



a 



(29) 
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where 



1 



and 



2 dA^dA, 
1 dx 



(30) 



2 dA^ 



(31) 



Ao 



The algorithm JT^ used to find the minimum utihzes a gradient search in the early 
stages of the fitting process, which transforms smoothly into a linearization of the fitting 
function as the fit converges. This is achieved by introducing a parameter A which sets the 
scale for the size of steps taken along the gradient. To do so, the diagonal elements of the 
curvature matrix a defined above are modified according to 

"/^/^ "m/. = (1 + ^) "m/^ ■ (32) 
Upon inspection of this equation, we note that if A 1, then 

SA^ ^ -^(3^ (33) 

and the search is approximately a gradient search. If, on the other hand, A <^ 1, then 

a'^^ ~ (34) 

and so the function is (approximately) linearized, as described in Eqs. ^ and |26|. The 
algorithm is designed such that the closer gets to a minimum (determined by the size of 
the changes in the fitting parameters required to 'step across' the minimum), the smaller 
the value of A chosen, and thus a transition from a gradient search to a linearization of 
is effected. 

At each scattering angle, a total of 10,000 randomly chosen sets of the A^'s were used 
as starting points to this algorithm. Approximately 95% of the time, the algorithm was 
successful and converged on a minimum in the multi-dimensional surface. It quickly 
became clear, however, that the function describing these observables contained many 
local minima, and the fit would often become 'trapped' in these shallower regions, rather 
than converging on the 'true' minimum value, Xbest- (This latter quantity was defined as 
simply the lowest value of Xmin determined in any of the 10,000 fits.) As a result, the 
algorithm was able to locate minima with Xmm — 2 xLst o^^Y about 10% of the time, an 
indication of the complexity of the space being probed. Nevertheless, this yielded a sample 
of roughly 1000 sets of amplitudes at each angle which gave reasonably good descriptions of 
all 16 measured observables. 

The procedures used to extract final values for each spin-operator amplitude will be 
presented in the next section. For now, we note that within these 1000 or so acceptable 
solutions, the magnitudes of the amplitudes, were generally found to be very stable, 

and largely independent of the choice of starting parameters. The values found for the 
corresponding phases, on the other hand, were highly dependent on the starting param- 
eters, and thus exhibited much larger fit-to-fit variations. Based on this, a second round 
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of fitting was undertaken, in wliicli the starting parameters were not chosen at random: in 
this analysis, the magnitudes were initially set to their best-fit values, as determined from 
the first round of fitting, while the phases were each stepped through all allowed values 
in a multi-dimensional grid search. Specifically, each phase was assigned a starting value 
between and Stt/S in steps of tt/3, with all possible combinations used as starting sets. 

As a final check on the robustness of the entire fitting procedure, we used a set of 
theoretical amplitudes (see Sect. VI) to generate values for each of the observables that had 
been measured in Ref. |T^, thereby producing a 'data set' comparable to that obtained in 
the actual measurement. Using this 'data' as input to our fitting code, we were indeed able 
to reproduce the input amplitudes, i.e., the algorithm could always find the correct solution. 
Once this had been established, we performed a more realistic test in which the calculated 
value of each observable was randomized with a root-mean-square deviation a equal to a 
typical experimental error bar for that observable. Using this data set, whose statistical 
precision matched that of the actual measurements, we were again able to reproduce our 
input amplitudes, within acceptable errors. 



V. RESULTS OF THE FITTING 

After performing 10,000 fits to the data at each scattering angle, our next step was to 
try to converge on a unique set of solutions for the Aj^'s. At each angle, we first discarded 
all fits in which the resulting minimum, Xmin' more than twice the lowest value 
found, Xbesf the number of degrees of freedom in our fitting function, this provided a 
confidence level of 85% that the true solution was included among the 1000 or so fits kept 



at each angle [0. Among these 1000 fits, though, the solutions tended to cluster tightly 
around a very small number of regions in parameter space, yielding a set of roughly 5-10 
distinct solutions for each angle. As mentioned previously, the amplitudes extracted from 
these different solutions were often quite similar in magnitude, but showed discrete (and 
correlated) variations in phase. 

To proceed further, i.e., to select from among these few distinct sets of solutions, it 
was necessary to introduce additional assumptions concerning the angle dependence of the 
amplitudes. By imposing the constraint that the amplitudes vary slowly and smoothly as 
a function of momentum transfer (as do all the observables described by these amplitudes), 
we were able to eliminate most of the remaining ambiguities in the values determined for 
the Aj^'s. We will illustrate this procedure with examples below. Before invoking such 
arguments, it is important to note that as the magnitude of a complex quantity goes through 
a minimum, both the real and the imaginary components must pass close to zero; hence the 
phase will typically change by roughly 180° as one passes through this minimum. Conversely, 
it is difficult to produce such a large phase shift unless a magnitude becomes very small. 

We now examine in detail the two amplitudes A^m and A^k- Shown in Fig. 1 are 
the fitted magnitudes and the relative phase of these two amplitudes, each plotted versus 
the center-of-mass scattering angle. (For plotting purposes, all phases lie in the range 
—180° < (pin < 180°.) Different plotting symbols at each angle correspond to different sets 
of solutions, and have been displaced slightly in angle. The magnitudes of both y4„„ and 
Akk are large over the range 5° < Oc.m. < 13° for all valid solutions, and decrease smoothly 
with angle. The phase difference (0„„ — (Pkk) is close to zero everywhere except for a single 
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solution at ^c.m. = 8.8°, indicated by a daggered 'X', which is near —180°. If one assumes 
that neither A„„ nor A^k passes near zero in this region, then this phase, and all other 
fitted values associated with this solution, are almost certainly unphysical. While one could 
attempt to make this argument more quantitative, e.g., by fitting the magnitudes of Ann 
and Akk with simple functions that did or did not pass near zero somewhere in this angle 
range, we feel that even a cursory examination of Fig. 1 tends to rule out the possibility of 
a zero-crossing. These arguments will become even stronger (though obviously somewhat 
more model-dependent) when we compare our results to a wide range of realistic theoretical 
predictions for these quantities, all of which display very smooth variations over this range 
of momentum transfer. 

It is interesting to note that the three quantities presented in Fig. 1 largely determine 
the observable Dog, which in the direction k • n = k • K = l/-\/2, k • q = has the value 

Doq '^{AnnA*j^j^) = |A„„||Axi^| sin(0„„ - (Pkk) ■ (35) 

Our measured value for this observable, being close to zero, drives the phase difference to 
either zero or 180°, but supplies little additional information to the fit. Constraints on 
the magnitudes of Ann and Akk, and the resolution of this 180° phase ambiguity, must 
therefore be supplied by other observables, demonstrating again the advantages of a large, 
diverse data set. 

By applying similar arguments to other Aj^'s, we were able to eliminate most of the 
remaining solution sets, and arrive at a nearly ambiguity-free determination of the magnitude 
and phase of each amphtude at all angles. For most quantities, the few distinct solutions 
left were sufficiently close in value that a simple average could be taken, with errors enlarged 
slightly to reflect the range of solutions included. Exceptions to this behavior occurred only 
at ^c.m. = 16.5°, where two of the phase differences exhibited two-fold discrete ambiguities 
that could not be resolved; an interpretation for this will be given in the following section. 
In all cases, the final values determined for the magnitudes and phases were consistent 
with those found in the fits which yielded the lowest minimum for x^, denoted here by x^est- 

As the final step in our analysis, the vahics obtained for Xbest were normalized to the 
number of degrees of freedom in the fit, to yield xl, a quantity which statistically should lie 
close to unity. In this work, the values for xl at 6*0.111. = 5.5°, 8.8°, 12.1°, and 16.5° were 0.75, 
2.14, 2.17, and 6.47, respectively. Because the fitting function should provide an accurate 
description of the data (that is, one is not gauging the appropriateness of the model in this 
case), a minimum xl in substantial excess of 1 suggests an underestimate of the input error 
for at least part of the fitted data set. To compensate for this, the error determined in the 
fitting process for the magnitudes and phases of each of the Aj^'s was artificially increased 

by a factor of \Jx^- This ensures that when the extracted amplitudes are used to determine 
best-fit values (with errors) for the observables, one will reproduce the measured input data 
within one standard deviation, on average. Thus, we believe that the errors quoted here for 
the magnitudes and phases of the A^i^^s provide a realistic estimate of the true uncertainties 
inherent in the fits to these data. 



13 



VI. COMPARISON OF AMPLITUDES WITH CALCULATIONS 



It has been shown previously |T^ that the large number of spin observables measured for 
this transition provides a severe test for any theoretical model. None of the calculations pre- 
sented in Ref. [jl2| could describe the momentum transfer dependence of all the observables 
over the entire range covered by the data. To determine more precisely where weaknesses 
lie in these models, it is useful to bring experiment 'closer' to theory by comparing not the 
measured observables, as was done in Ref. [|12| but by comparing theoretical predictions for 
the Aifj_ amplitudes with values extracted directly from the data, as described here. 

In this work, we will consider five different sets of predictions for the scattering ampli- 
tude. To carry out these calculations, many details must be specified, such as: the method 
used to generate the distorting potential, and how well this describes the elastic scattering 
data; assumptions made regarding the structure of the excited state; the extent to which 
medium corrections are incorporated into the effective interaction; the handling of exchange 
contributions to the interaction; etc. Here, we will examine predictions for the Aj^'s derived 



from the same five models as those to which the observables were compared in Ref. |]12 
Extensive discussion on the content of each calculation is presented in that paper, and the 
interested reader is encouraged to refer to Ref. for more detail; only fairly brief descrip- 
tions will be provided here. 

Of the five calculations to be shown, all use free t-matrices for the effective interaction, 
and Cohen-Kurath matrix elements to specify the transition to the excited state. Two 
are relativistic calculations in which both direct and exchange contributions are included 
explicitly (DREX) The distorting potential for the incident and outgoing projectile 

waves, however, is generated either "self-consistently" (the same interaction that induces 
the 0+^1^ transition is also folded with the ^^C ground state transition density to pro- 
duce the distortions), or from an optical potential, with parameters fit to p+^^C elastic 
scattering cross section and analyzing power data. In all of the figures, the predictions 
of these two calculations will be shown as a solid line (self-consistent) and a long-dashed 
line (optical potential). Two nonrelativistic calculations also include both direct and ex- 
change contributions explicitly (DW81) |TB|, and also incorporate the effects of distortion 
using the same two methods, i.e., either self-consistently (dotted line) or from the same 
optical potential as was used in the relativistic calculation (dot-dashed line). Finally, we 
will compare our values of the scattering amplitudes to those predicted using a relativistic 
calculation in which the full interaction has been parameterized in terms of direct scat- 
tering processes only, with distorted waves generated self-consistently (short-dashed line). 
Although exchange processes are expected to contribute significantly to this reaction, this 
last calculation appeared to describe the observables measured for this reaction better than 
any model in which exchange contributions had been explicitly included |jl2|. This would 
suggest that most current methods of accounting for exchange processes are inadequate to 
describe proton-nucleus scattering at intermediate energies, at least for unnatural parity 
transitions. 

Before making detailed comparisons between data and theory, there is one ambiguity 
which must be pointed out. Because the individual spin-operator amplitudes have dimen- 
sions of /ib/sr^, cross section data |T3| were used to provide an overall normalization factor 
for all spin observables considered in the fitting procedure. This ensures, for example, that 
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Eq. H is obeyed, and the sum of the squared amphtudes equals the measured cross section. 
However, as can be seen in Fig. 2, the five calculations considered here do not reproduce 
this cross section. In particular, all of the calculations underpredict da/dflp at small angles, 
except for the DW81 calculation using distortions generated in an optical model, which 
overpredicts the data at small angles, then decreases much more rapidly with angle than the 
data. The DW81 calculation which uses self-consistent distortions, on the other hand, is a 
much more shallow function of angle than the data suggest. 

In light of these discrepancies, one could consider scaling all of the calculated magnitudes, 
l^j^l, by a single multiplicative factor, in order to more closely reproduce the measured 
cross section. In this way, one is effectively comparing the relative sizes of the spin-operator 
amplitudes for each model to the measured values. This also ensures that the dominant 
amplitudes will be well reproduced. On the other hand, the differences observed between 
the predicted and measured cross sections might result from particular amplitudes (especially 
the larger ones) being grossly over- or under-predicted, while the calculated values for others 
are actually in close agreement with the experimentally determined values. In this case, one 
is best served by not applying an artificial normalization, and directly comparing data 
and theory for each amplitude. Because our primary goal is to identify more narrowly 
the weaknesses in individual calculations, we have adopted the latter approach here. The 
potential drawback, of course, is that a calculation that is correct in all respects other 
than reproducing the measured cross section will systematically miss each of the extracted 
amplitudes by roughly the same factor. 

With the above caveat in mind, we now compare the predictions of these five models 
for the magnitude and phase of each amplitude to the values deduced from our fits. The 
quantity Ano, shown in Fig. 3, is the amplitude for polarizing the recoil ^^C nucleus along 
the n direction (perpendicular to the scattering plane) with an unpolarized incident proton. 
Its magnitude (upper plot) is predicted by each calculation to be quite weak over the entire 
angle range studied in this work, a feature also seen in the data. The differences among the 
various calculations are, in absolute terms, very small. Thus, given the level of precision with 
which this quantity can be determined experimentally, not much useful information can be 
obtained from |v4„o|, other than confirming the small probability for producing this particular 
spin configuration. Due to the small size of \Ano\, the phase difference 0„o ~ (pnn, shown 
in the lower half of Fig. 3, exhibits a twofold discrete ambiguity at the largest scattering 
angle. As discussed earlier, if the magnitude of a complex amplitude passes near zero, its 
phase can change by nearly 180°. Our data suggest that A^o passes near zero (has a local 
minimum) somewhere around ^c.m. = 16.5°, but our measurements can not establish this 
unambiguously. We also note that at smaller angles the phase difference is relatively fiat, 
in agreement with all of the calculations, although the values deduced from the data are 
significantly more negative (by ~ 90°) than any of the calculations predict. 

In contrast to this weak amplitude, we next examine Ann, shown in Fig. 4, which is the 
amplitude for polarizing the recoil nucleus along the n direction when the incident proton 
is also polarized along fi. The data show that is large throughout the angular range 

studied, and decreases smoothly with increasing momentum transfer. All five calculations 
predict this general behavior, but differ significantly in strength, relative to the precision 
with which this quantity has been determined. It is useful to note the striking similarities 
between this figure and Fig. 2, the unpolarized cross section. (For a more quantitative com- 
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parison, one would need to square the results shown in Fig. 4.) Just as for da/dflp, the 
values of at small momentum transfer are underpredicted slightly by four of the cal- 
culations, and overpredicted by the DW81 calculation which uses optical model distortions. 
The momentum transfer dependence of the data is reasonably well described by the three 
relativistic calculations, but is too steep (shallow) for the DW81 calculations that use optical 
model (self-consistent) distortions. The inabilities of the models to reproduce this particular 
amplitude are thus directly reflected in the discrepancies found between the predictions and 
measured values for the scattering cross section. 

Similar behavior is seen in Fig. 5 for the other large amplitude, Akki the amplitude for 
polarizing the recoil nucleus along the average momentum direction using an incident proton 
polarized along the same direction. The measured magnitude l^xi^l (upper plot), like 
closely follows the angular dependence of the differential cross section, and decreases mono- 
tonically with 6'c.m.- In this case, the three relativistic calculations quantitatively describe 
this behavior, and agree with the data at all values of momentum transfer. The two non- 
relativistic calculations, on the other hand, do very poorly, either over- or underpredicting 
the strength of this amplitude at small angles, and predicting angle dependences which are 
either too steep or too shallow, depending on the method used to generate distortions. Also 
shown in Fig. 5 is the phase difference 0„„ — (I)kk^ which is predicted by all of the calcu- 
lations to be very close to zero over this entire angular range. The data support this idea, 
albeit with a fairly large statistical uncertainty, suggesting that the two largest amplitudes 
contributing to this reaction are closely matched in phase. 

The third "diagonal" amplitude, Aqq, has some intriguing properties. Aqq is the ampli- 
tude for polarizing the recoil ^^C nucleus along the q (momentum transfer) direction, with 
the proton also polarized along q, i.e., this amplitude is associated with the spin-operator 



combination (S-q)((T-q). In a meson-exchange formulation of A^A^ scattering |TP[, the one- 
pion exchange amplitude is associated with the spin-operator (o"i-q)((T2-q), which is very 
similar. Perhaps of more significance, in a relativistic PWIA, Aqq is the only amplitude that 
contains contributions from the pseudovector invariant A^A^ amplitude, which carries the 
same quantum numbers as the pion. Because the excitation of the 15.11 MeV state is a 
AT = 1 transition, one-vr and one-p exchange should dominate the reaction process. The 
potentials for vr and p exchange interfere constructively at very small momentum transfer, 
but the one-vr exchange term changes sign relative to one-p exchange at relatively low mo- 
mentum transfer ||2^ , which can lead to large cancellations. Thus, a zero crossing in the Aqq 



amplitude may be a direct reflection of vr — p interference in the reaction mechanism for this 
transition. 

Inspection of the extracted values for the magnitude and phase of Aqq, plotted in Fig. 6, 
clearly indicates that this amplitude does indeed cross zero somewhere between 6'c.m. = 
5.5° and 12.1°. The magnitude \Aqq\, which exhibits a very different q dependence than 
either \Ann\ or is very well described by all of the relativistic calculations, which 

each predict a zero crossing somewhere near 6'c.m. ~ 9°. The data are not described quite 
as well by the nonrelativistic calculation which uses self-consistently generated distortions, 
and are in strong disagreement with the other DW81 calculation. The large experimental 
uncertainty in the phase difference at 6'c.m. = 8.8° supports the idea that a zero crossing 
occurs near this point. The two relativistic calculations which have distortions generated 
self-consistently describe the momentum transfer dependence of both the magnitude and 
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phase of Aqq quite well, and thus best describe the physics contained in this amplitude. 

The last two amplitudes, the off-diagonal terms A^q and Agx, are also intriguing. Physi- 
cally, they represent the amplitudes for polarizing the recoil ^^C nucleus along either the K or 
q direction when the proton probe is polarized along either the q or K direction, respectively. 
In a nonrelativistic PWIA, these two amplitudes should be identically zero for this transi- 
tion, and only become non-zero if non-local effects, such as knock-on exchange, are explicitly 
included |^ . In a relativistic formulation, on the other hand, Akq and Agx include contribu- 
tions due to linear couplings between the upper and lower components of the bound nucleon 
wave function, even in PWIA. Because the lower components are momentum-dependent, the 
nucleon is manifestly off-shell, and non-local effects arise "naturally." Formally, A^q and 
AqK are proportional to the tensor and axial vector pieces of the invariant NN amplitude, 
respectively. By carrying out the spin algebra 0, these can be written in terms of the 
composite spin-convection current amplitudes {a x J) and {a ■ 3), again respectively. Thus, 
these two amplitudes should be sensitive to the off-shell behavior of the nucleons inside the 
^^C nucleus. 

In Fig. 7, we see that the extracted magnitude of Akq is significantly smaller than all 
of the predicted values, especially at small angles. This suggests a quenching of the tensor 
component of the A^A^ amplitude within the nucleus, though one can not tell if this is 
a nuclear structure or an interaction effect. It is very curious that the only relativistic 
calculation that does not include exchange (short-dashed line) also predicts the largest 
magnitude for an amplitude that, at least nonrelativistically, is driven largely by exchange 
contributions. Unfortunately, little insight is provided by the phase difference, 0nn — 4>Kq- 
Note, however, that if the phase found at the largest angle is simply shifted by 360°, the 
general shape of the angular dependence of the phase difference is followed reasonably well 
by three of the calculations, albeit with a ~ 90° offset. 

In striking contrast, the other off-diagonal amplitude, Agx, is predicted to be somewhat 
smaller than Akq over most of this angle range, while the experimentally determined values 
are seen in Fig. 8 to be quite a bit larger. Only the DWIA calculation using optical model 
distortions is consistent with the typical strength suggested by the data. Figure 8 also shows, 
though, that the limited statistical precision with which this amplitude has been determined 
makes it difficult to draw any conclusions about the momentum transfer dependence of this 
quantity. This limitation is also evident in the phase difference, in which a tight clustering 
about zero is broken only at the largest angle, where a two-fold discrete ambiguity is observed 
at 16.5°. As was the case for A^o, our statistical accuracy is such that we cannot determine 
experimentally whether this amplitude is crossing zero in this angle regime or not. Two 
calculations predict amplitudes that, in the complex plane, pass near zero on one side (A0 = 
-1-180°), two pass on the other side (—180°), and one phase remains constant. All of these 
possibilities are consistent with one of the allowed solutions extracted from the data. 

VII. SUMMARY AND CONCLUSIONS 

We have made a model-independent determination of the full transition amplitude for 
the ^'^C{p,p'y^C* (15.11 MeV, 1+, T = 1) reaction at an incident beam energy of 200 MeV at 
four scattering angles. This represents the first such determination for any hadron-induced 
nuclear transition other than those with Ji = Jf = 0. By imposing only loose "smoothness" 
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constraints on the momentum transfer dependence of the individual spin-operator ampli- 
tudes, we have performed a nearly ambiguity-free extraction of these quantities, which has 
provided deeper insight into the physics driving this transition. 

As expected theoretically, the two diagonal spin-operator amplitudes A„n and A^k were 
found to be the dominant amplitudes over the entire momentum transfer range studied. Each 
of these was better described by the three calculations carried out in relativistic frameworks 
than by the two nonrelativistic calculations we considered. The third diagonal amplitude, 
Agg, which has a spin-operator structure similar to that of one-pion exchange, exhibited 
behavior characteristic of a "zero crossing." The value of momentum transfer at which this 
occurred was again much better matched by the relativistic calculations than the DWIA 
predictions. This crossing has physical significance in that it may reflect interference between 
TT and p exchange, and hence may serve as a gauge of the relative strengths of these two 
contributions within the nuclear medium. 

The three off-diagonal amphtudes are all much weaker than the three just discussed, 
and as such were determined with much larger experimental uncertainties. The magnitude 
of the Ano amplitude provided little insight, though the phase was consistently (i.e., at 
all angles) off by about 90° relative to most of the predicted values. The amplitudes Axq 
and AgK are expected to be sensitive to the off-shell behavior of the nucleons inside the 
^^C nucleus, and should therefore probe the non-local or exchange nature of the scattering 
process. Despite the sizable errors on the experimentally determined values for these two 
amplitudes, neither are described well by any of the five calculations considered here, which 
may indicate problems in our present treatment of non-local effects in both relativistic and 
non-relativistic frameworks. 
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FIGURES 



FIG. 1. Sets of fitted values for the magnitudes of Ann (top) and Akk (middle), and their 
relative phase difference (bottom), each plotted vs. ^c.m.* I^ each case, different plotting symbols 
at each angle represent the results of different solutions. The symbols have been displaced slightly 
in angle for (some) clarity. 

FIG. 2. Differential cross section for the ^^C{p,p')'^^C* (15.11 MeV) reaction at 200 MeV. The 
data are from Ref. [13]. The five curves shown are described in the text. 

FIG. 3. Magnitude (upper) and relative phase (lower) of the amplitude Ano plotted vs. the 
center of mass scattering angle. The phase difference is with respect to the Ann amplitude. The 
error bands represent best fit solutions to the data of Ref. [12] . The five curves shown are theoretical 
predictions described in the text. 

FIG. 4. Same as Fig. 3, but for the amplitude Ann- In our fitting procedures, the phase of this 
amplitude was defined to be zero. 

FIG. 5. Same as Fig. 3, but for the amplitude Akk- 

FIG. 6. Same as Fig. 3, but for the amplitude Aqq. 

FIG. 7. Same as Fig. 3, but for the amplitude Akq- 

FIG. 8. Same as Fig. 3, but for the amplitude AqK- 
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